Universal rates for reactive ultracold polar molecules in reduced dimensions 
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Analytic expressions describe universal elastic and reactive rates of quasi-two-dimensional and 
quasi-one-dimensional collisions of highly reactive ultracold molecules interacting by a van der Waals 
potential. Exact and approximate calculations for the example species of KRb show that stability 
and evaporative cooling can be realized for spin-polarized fermions at moderate dipole and trapping 
strength, whereas bosons or unlike fermions require significantly higher dipole or trapping strengths. 
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The successful production of relatively dense gases [T] 
or lattices of ultracold molecules in their ro- vibrational 
ground states opens up a number of new opportunities 
in physics and chemistry [5J 0]. Reaction rates of ul- 
tracold molecules can be quite large, as measured and 
calculated for the fermionic species 40 K 87 Rb [5--7J. This 
highly reactive species belongs to the class of molecules 
that have universal reaction rates that can be calculated 
analytically from a knowledge of the long-range van der 
Waals (vdW) potential alone, given unit probability of 
reaction at short range [5J . When an electric field is used 
to polarize a gas of this polar species, reaction rates be- 
come even larger and significantly limit the lifetime of 
the dipolar gas 18] . A necessary condition for achieving 
stable dense samples of such reactive ultracold molecules 
is to be able to understand and control such loss pro- 
cesses. This is essential, for example, to achieve evapo- 
rative cooling in order to reach quantum degeneracy and 
thus utilize the molecules for new and exotic condensed 
matter phenomena that have been proposed [9HTT] . Ref- 
erence Q~2] predicted that molecular dipoles can be stabi- 
lized by tightly confining the molecules in a single plane 
of quasi-two-dimensional (quasi-2D) geometry. This is 
because the molecules can be arranged so as to experi- 
ence repulsive dipolar forces at long range and thus never 
come together within chemical interaction distances in 
order to react. First steps have been taken towards un- 
derstanding these effects by recent quasi-2D calculations 
with adiabatic potentials 1 13] or quantum dynamics for 
fermionic 40 K 87 Rb [14]. 

Here we characterize quasi-2D and quasi-lD elastic and 
reactive collision rates of a broad class of bosonic and 
fermionic molecules that have universal rate constants. 
Analytic expressions apply in the case of vdW interac- 
tions. Quantum dynamical and approximate calculations 
for quasi-2D collisions of polar molecules show how colli- 
sion rates scale with confinement length, dipole strength, 
and collision energy. Our analysis connects the critically 
important experimental domain between the vdW and 
strong dipole limits, and provides a general framework 
for future research on molecular cooling, chemical reac- 
tivity, and condensed matter applications for species with 



universal collision rates. 

Universal collision rates occur for chemical species that 
have near unit probability of reaction or inelastic relax- 
ation when they are close enough together, on the order of 
typical chemical interaction distances, a c < 1 nm, where 
strong chemical forces permit the reaction to occur [BJ. 
Universal collision rates arc completely determined by 
quantum threshold dynamics associated with the long 
range potential. Consequently, there are two distinct 
classes of mixed alkali-metal diatomic species. Those 
with energetically allowed reaction channels, KRb, LiNa, 
LiK, LiRb, and LiCs jTHj, are expected to be universal. 
By contrast, the species NaK, NaRb, NaCs, KCs, and 
RbCs have no reactive channels [15] , and are expected to 
be non-universal in their ground rotational, vibrational, 
and spin state. However, even these species, as well as 
nonreactive homonuclear dimers, can have universal in- 
elastic relaxation rates when vibrationally or rotationally 
excited [TH [17] . Universal species do not have scatter- 
ing resonances, since threshold bound states decay too 
fast to make one, whereas non-universal species can have 
numerous scattering resonances. This paper treats the 
broad class of universal collisions. 

Figure [l] shows the long-range potential V(p, z) for two 
dipoles in quasi-2D geometry and illustrates essential fea- 
tures of reduced dimensional collisions, where 
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Here r = (p, </>, z) represents the distance between the 
two molecules in cylindrical coordinates, and r = |r|. 
The molecules are confined in the z-direction by a har- 
monic trap of frequency O and characteristic length 
^ = y/h/pCl, where p is the reduced mass of the 
pair. The dipoles are assumed to be aligned along z, 
so the projection m of their relative angular momen- 
tum is conserved. The second term represents the m- 
dependent centrifugal potential in 2D. The third term 
is the isotropic vdW potential, assuming the molecules 
are in their rotational ground state, with a vdW length 
a = [27r/r(l/4) 2 ] (2^C 6 /h 2 ) 1/4 QUE]. The last term 
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FIG. 1: Potential V(p, z)/h£l versus p/a^ and z/ah for m = 
1 for 40 K 87 Rb in a {1/2% = 50 kHz trap, where a h = 56.4 nm, 
a = 6.25 nm with Cq = 16130 a.u. 6, 19 . (a) vdW only case, 
d — 0, (b) dipolar case for d = 0.2 D (with ID = 1 Debye = 
3.336 10 - 30 Cm). Saddle points (white circles), (a) in-plane 
and (b) out-of-plane minimum action paths (dashed lines). 
Dotted half-circles indicate a and a d . 



is the anisotropic dipolar potential, with induced dipole 
moment d and dipolar length a<j = [id 2 /h 2 . 

Our model relies on the separation of length scales: 
o-k 3> &h 3> d 3> a c , where a K — 2tt/k is the DcBroglic 
wave length for a collision with relative kinetic energy 
E K = h 2 n 2 /2[i. These inequalities are readily satisfied 
for experiments with 40 K 87 Rb , for which a K , and 
a can be on the order of hundreds of nm, tens of nm, 
and less than 10 nm, respectively. Under these condi- 
tions, collisions are essentially quasi-2D HI], or if 
additional confinement at frequency f2 is provided along 
x, quasi-lD [22]. Figure [T|a) illustrates the quasi-2D 
vdW case with d = for |m| = 1, where the minimal- 
action path followed by the two colliding particles lies 
in the plane z = of the contour diagram of V(p,z). 
The competition between the centrifugal barrier and the 
vdW attraction determines a single saddle point at a dis- 
tance r > a, separating the long-distance 2D-scattering 
region from the short-range "core," r < a, where the 
molecules are accelerated towards one another by the 
attractive potential and experience 3D scattering. Two 
additional out-of-plane saddle points appear when d in- 
creases so a d /a > 2.71(a/dh) 3 / 2 , as in Fig. [ijb), indicat- 
ing a crossover to dipolar-dominated scattering, which is 
fully reached for ad > ah S> d. We show below that a 
strong enhancement of the ratio of elastic to reactive col- 
lisions in this regime will allow for an efficient cooling of 
the molecular gas. 

Reference [5J characterized elastic and reactive col- 
lision rates for 3D collisions by a complex scattering 
length dj{k), where hk is the 3D momentum. For the 
special universal class of highly reactive molecules with 
unit short-range reaction probability, ao(fe) = (1 — i)a 
for s-wave collisions of like bosons or unlike fermions, 
and a\(k) = (—1 — i){kd) 2 d\ for p-wave collisions of like 
fermions, with ai w 1.064a [6 . These explain the mea- 



TABLE I: Lengths Lj(n) and ratios rj(n) defining the 
universal complex scattering phase in N dimensions; £ is the 
Riemann C-function, B = 0.905 HJ and W(0) = 0.328 [21]. 
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sured rates of 3D collisions of ultracold 40 K 87 Rb [3 [6] . 
Universal species have only incoming scattering current 
in the entrance channel in the vdW "core" of the collision, 
r < a, and this provides a universal boundary condition 
for both the vdW and dipolar cases illustrated in Fig. [l] 
While the general model can be extended to nonreactive 
species like RbCs, where one expects a rich resonance 
structure, we here focus on universal species like KRb. 

Applying universal theory to quasi-lD and quasi-2D 
collisions in a vdW potential is straightforward for the 
d = case by combining the methods and notation of 
Refs. [6l|23]. We give here only the resultant formulas, 
valid for nd <C 1. Assume the molecule is prepared in its 
vibrational, rotational, and spin ground state and in the 
ground state of confined motion. Only the first channel j 
of a coupled channels expansion is needed near threshold 
for small d, and we set the index j — for like bosons 
or unlike fermions and j = 1 for like fermions. Elastic 
and reactive collisions in N dimensions, N — 1,2,3, are 
described by an S'-matrix element Sjj — exp(z^) written 
in terms of a complex phase 0j(k) with 
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where n = p,q,k represents the momentum in ID, 2D, 
and 3D respectively 23 . The quantity dj(n) on the right 
hand side defines the complex scattering phase through 
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Tabled gives the lengths Lj(n) and ratios rj(n) for N — 
1,2,3, where £ = d/a^ and £i = did 2 /a\. For the case 
in Fig. [lj Co = 0.111, £i = 0.00145, q = (1/126) nm- 1 , 
q 2 a\ = 0.200 for E K = A: B 240nK, where k B is Boltz- 
mann's constant. Assuming d/a^ <C 1, we find rj <C 1 
and the right-hand factor in Eq. ^ can be approximated 
by (— iy — i, as in 3D [B]. The lengths Lj change only 
weakly across dimensions N, so the scaling is mainly 
given by the factor a^~ N . 

The elastic JCf and reactive K.™ scattering rate con- 
stants follow from the formulas in Ref. 231: 
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where go = 1/tt, 2/tt, 2 and gy = 1/tt, 4/7r, 6 for molecules 
colliding in like spin states in N = 1, 2,3 dimensions re- 
spectively, which take into account the iV-fold degeneracy 
in j = 1. Unlike bosons or fcrmions have rate constants 
(/Co + /Ci)/2. The upper bound, or unitarity limit, on the 
rate constants follows immediately upon replacing the S- 
matrix expressions by their upper bounds, |1 — 5^1 2 < 4 
and 1 — |SV,-| 2 < 1, respectively. The elastic Tf = tCfn 
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(b) reactive (j=l) 
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with n the density of the collision partner in N di- 
mensions (units of cm~ N ), n/a^~ N an equivalent 3D 
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AT = 1 or 2 dimensions 
-Ar gives the ratio of elas- 



FIG. 2: Quasi-2D rate constants per energy, K\/(E q /k&), for 
elastic (a) and reactive (b) collisions of identical KRb fermions 
(j — 1) versus dipole strength d for an Q/2n = 50 kHz trap 
(cf. aa = ah for d = 0.25 D) and different energies Eg/ky^; CC 
(solid lines) and UBA/instanton (dashed lines). 



density (units of cm J ) for 
and rij(K) = 2Lj(k)k n ^ 2 / 'a^' 

tic to reactive collision rates. The factor fj(n) — 

[l + r,j (k) + r.j (k) 2 /2 + rjj (k) + rjj (k) 2 /2] 1 approaches 
unity as k — > 0, except for j = for N = 1,2, where 

fo(p) -> [pa h /^o] 2 /8 and f (q) -> 7r/21n 2 \2B/ir?a£[, probability P? D for tunneling through the barrier sepa- 
respectively. However, for realistic traps and energies in 
the nK-regime the full expression for fo(q) is required in 
2D, since the logarithmic term only becomes dominant 



rating the long-distance 2D-scattering region r>a from 
the short-range vdW core r < a (see Fig. [lj : 



at much lower energies. The expressions in Eqs. (5a) and 



(5b) give explicitly the scaling with k and ah, as well as 
all the known threshold law limits for collision rates in 
reduced dimensions [2T)H2"2l |2"S"] . 

We use both coupled-channels (CC) methods and an- 
alytic/semiclassical approximations to show the effect of 
the dipole moment on 2D collision rates. The former uses 
a spherical harmonic basis set and a renormalized Nu- 
merov method to propagate the wave function ^(r) with 
universal incoming wave boundary conditions in the vdW 
core [6] out to distances r ^ ah. Then ^(r) is matched 
onto a cylindrical basis and propagated to larger p to 
yield the 2D S'-matrix. The numerical results agree with 
the vdW limits from Eqs. Q, and Tableland the 
dipolar results in Ref. [T4"] . 

Elastic collisions are well-described by a unitarized 
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Born approximation (UBA), S 
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iK 2 ?) 1 , where the K- matrix element in 2D includes 



jj 



the vdW term from Eq. ^ plus the dipolar term, 

" 2D (g) = -5f (?) + 2^Mqa h ) , (6) 



with (f> {x) = -0.65471 + 0.94146a; - 0.39010a; 2 + 0(x 3 ) 
and <j>i(x) = -0.35555x + 0.36042x 2 -0.13417x 3 +O(:r 4 ). 
Equation ^ shows how the elastic rate constant scales 
with q, ad, and ah- Figures ^ a) and|3ja) show that the 
UBA is an excellent approximation for dipolar bosons 
and fermions for a wide range of realistic E q , fl and d. 

Reaction rates can be estimated using an instanton 
technique [T21 HE] with V(p, z) to get the transmission 



Here, r^ 1 is the path of minimal action given by the classi- 
cal trajectory of a particle in the inverted potential with 
inner and outer turning points, r x and r 2 , and Sf the 

associated Euclidian action. We use Vj = Vj + h 2 /8pp 2 
to take into account the semiclassical Langer-correction 
to the centrifugal term in the potential (m 2 — 1/4 — » m 2 ), 
insuring correct threshold laws for P 2D (q). We find 
Ay = 1.0297(a/ah) by equating the analytic expression 
for P 2D (q) to the analytic 1 - l^nl 2 in Eq. Q in the 
vdW limit d — > 0, assuming Aj is independent of d for 
a d < Oh- Since Sq — > as d — > due to the disappear- 
ance of a centrifugal barrier for j = 0, we set Ao = 1 to 
ensure unitarity is satisfied and only use Sq to estimate 
/Cq 1 at finite d where a barrier exists and ad < ah- 

Figure [2] shows that the variation of JCf /(Eg/hs) and 
K. T i /(Eq/ks) with energy is relatively weak. As in 3D, 
IC T i /(Eq/ks) is independent of E q at low d, but quite 
unlike in 3D [7j |8] , it decreases relative to the vdW limit 
when d increases. The instanton method gives the qual- 
itative explanation. At small d the barrier to the in- 
plane path increases with d, thus decreasing K.\ e . As d 
increases, the existence of out-of-plane saddle points gives 
alternative paths with a lower barrier, so JC\ C starts to in- 
crease. As d increases more, the increasing out-of-plane 
barrier strength eventually will cause P 2D (q) to again 
decrease with increasing dipole, evident in the CC and 
instanton 150 kHz strong trap case in Fig. Ub). Com- 
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FIG. 3: Quasi-2D elastic (el) and reaction (re) rate constants 
Kj for identical KRb (a) bosons ( j = 0) and (b) fermions (J = 
1) at a collision energy of E q = 240 nK = h 5 kHz for three 
different trap frequencies Q/2tv = 25 kHz, 50 kHz, 150 kHz. 
CC (thick lines), UBA/instanton (thin lines), unitarity limits 
(dotted lines). Vertical arrows indicate where rjj = 100. 



kHz trap for d = 0.3 D. Thus, stability and evaporation 
for universal polar bosons may be achievable. 

In conclusion, we have developed universal analytic ex- 
pressions for quasi-2D or quasi-lD collisions of highly 
reactive ultracold bosonic or fermionic molecules in the 
absence of an electric field and have calculated quasi-2D 
collisions of universal polar molecules with a dipole mo- 
ment. While prospects for stability and evaporative cool- 
ing in 2D are much better for universal polar fermions 
than bosons, either species would benefit from larger 
dipole moments or tighter confinement. Non-reactive, 
non-universal species are expected to be very different 
from universal ones and to show numerous shape or Fes- 
hbach resonances as electric or magnetic fields are tuned. 
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and the EU NAME-QUAM. We thank G. Quemener, J. 
L. Bohn, and Jun Ye for discussions. 



paring Figure |3ja) and[3jb) shows the instanton approx- 
imation gives the qualitative trends for the boson case 
even better than for the fermion case at finite d. 

Figure [3] illustrates the stability and cooling proper- 
ties expected for universal polar bosons and fermions. In 
both cases /C^ 1 changes from the vdW limit by rapidly in- 
creasing with d until it approaches the 2D unitarity limit 
at large d. Stability requires that T¥*(q) remain small 
enough that the lifetime 1 /T™ (q) is sufficiently long, or- 
der of 1 s or longer, as achieved for 40 K 87 Rb fermions 
in 3D [SJ. Equation (5b I predicts that the reaction rate 



per particle for identical fermions in the vdW limit is only 
times lower in 2D than 3D, if the 2D system has the 
same equivalent 3D density n 2D /ah- The initial decrease 
in IC\ C indicates not only increased fermionic stability in 
2D, but also the possibility of evaporative cooling, which 
requires the ratio rjj = Kf/Kf ^> 1. Figure pi b) shows 
that 7/i reaches a magnitude near 100 for a — 0.18 D, 
nearly independent of the trap O. Figure [2] shows that 
the ?7i ratio increases with lower E q , indicating that the 
evaporation improves as the 2D gas cools. 

Figure [3^ a) shows that there is less room to improve 
the bosonic elastic rg^g) with increasing d, since it is only 
an order of magnitude below unitarity in the vdW limit. 
The reactive Tq ^) at low dipole strength is much larger 
than for fermions, so universal polar bosons have shorter 
lifetimes and poorer stability than fermions at the same 
dipole and trap strength. Furthermore, getting ryo > 100 
requires either large d or large f2. Figure [3] shows that 
r/o = 100 in a 50 kHz trap near d = 0.5 D or in a 150 
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